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ABSTRACT 


Nonlinearities  are  very  significant  for  the  surging 
motion  of  a  ship  in  a  following  sea.  Results  are  presently 
available  only  for  the  surging  motion  in  regular  following 
waves.  The  law  of  superposition  cannot  be  used  to  treat  the 
nonlinear  surging  motion  in  an  irregular  sea.  Therefore, 
it  is  not  possible  to  extend  the  results  obtained  in  regular 
waves  to  predict  the  behavior  in  an  irregular  sea. 

A  method  to  overcome  these  difficulties  to  invest¬ 
igate  nonlinear  surging  motion  in  an  irregular  sea  is  pre¬ 
sented  in  Part  1  of  this  report.  Using  this  method,  results 
are  obtained  for  nonlinear  surging  motion;  in  particular, 
solutions  are  given  for  the  acceleration  of  a  ship  to  an 
increased  speed  and  for  the  subsequent  run  at  that  increased 
speed. 

Based  on  the  results  obtained  in  Part  1,  the  problem 
of  broaching  is  discussed  in  Part  2.  A  connection  between 
broaching  and  nonlinear  surging  motion  was  presumed  in  previous 
discussions  about  broaching.  This  connection  is  confirmed 
and,  furthermore,  nonlinear  surging  motion  is  found  to  be 
the  most  important  pre-condition  for  broaching. 
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SURGING  MOTION  OF  A  SHIP 
IN  AN  IRREGULAR  FOLLOWING  SEA 


INTRODUCTION 

To  compute  the  surging  motion  of  a  ship  in  a  head 
sea,  a  linearized  equation  of  motion  can  be  assumed.  This 
leads  to  a  straightforward  solution  in  regular  waves.  The 
surging  motion  in  an  irregular  head  sea  can  be  treated  as  a 
random  process  and  solved  by  the  method  of  St.  Denis  and 
Pierson. 1 

The  assumption  of  linearization  may  not,  however, 
be  valid  for  surging  motion  in  a  following  sea,  since  non- 
linearities  become  significant  in  this  case.  This  fact  wa3 
demonstrated  in  refs.  2  and  3  where  it  was  noted  that  in  a 
certain  range  of  propeller  thrust,  the  oscillatory  surging 
motion  vanishes  in  regular  following  seas  and  the  ship  is 
forced  to  run  at  the  same  speed  as  the  wave.  This  behavior 
is  very  striking  and  is  believed  to  be  connected  in  some  way 
with  many  known  cases  of  ships  broaching  in  following  sea3. 

Because  this  striking  behavior  in  a  regular  follow¬ 
ing  wave  is  due  to  nonlinearities,  no  conclusion  about  the 
analogous  behavior  in  an  irregular  following  sea  is  possible. 
A  treatment  of  the  surging  motion  in  an  irregular  following 
sea  must  include  the  influence  of  these  nonlinearities  and 
great  difficulties  arise  from  this  requirement.  No  scien¬ 
tific  treatment  of  this  Bort  is  known  to  the  author. 

This  paper  attempts  to  overcome  these  difficulties. 
No  complete  solution  of  the  problem  is  attained,  but  some 
answers  are  found  —  particularly  to  the  question  of  whether 
the  ship  can  be  accelerated  to  a  greater  extent  than  may  be 
predicted  by  a  linear  equation  of  motion. 
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THE  EQUATION  OF  SURGING  MOTIONS 
AND  THE  ROLE  OF  NONLINEARITIES 


A  two-dimensional  irregular  sea  may  be  represented 


by: 


h(x,t)  =  to-f  r(<-on)Aton  cos  pin2  x  +cunt+un  j 


or  = 


J00  _ 

o  cos|^i!l-x+cot+  a(oc)j  ^r(o)) 


(1) 


dew 


The  force  in  the  longitudinal  direction  produced 
by  a  regular  longitudinal  wave  ( the  exciting  surging  force) 
can  be  approximated  by  using  the  Froude-Kryloff  hypothesis. 
The  ratio  of  this  force  divided  by  the  wave  amplitude  can  be 
represented  by: 


I  =27T.y.f  X,  =i£f  V.f(H^%  (2) 

h  x.  VL'  g  g  ' 

where  f  (£)  =  f(iifi-)  and  denotes  the  forcing  function  of  the 
g 

exciting  surging  force.  This  function  depends  on  the  form 
of  the  ship  and  on  the  ratio  of  wave  length  to  ship  length. 
For  the  numerical  computation  cai’ried  out  and  represented  in 
this  paper,  the  followihg  is  used  as  a  possible  function 
(Fig.  1): 

f(|)  -  3  sln  (gr°  ~  ^  cos 

/CUgL\  3  (3) 

(Sg"J 

In  the  same  manner,  the  surging  force  on  a  ship  in  a  two- 
dimensional  irregular  sea  running  longitudinally  becomes: 


00 
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It  is  convenient  to  divide  this  force  by  the  mass  of  the 
ships  m  =  ~*.  The  result  has  the  dimensions  of  acceleration: 

o 

-j|  =  )0  cosjcp  X  +  cut  +  a (cl>) J  ^ cu4f2(^)  r(cb)  dec  (5) 

If  we  assume  that  the  ship  is  moving  at  a  constant 
mean  velocity  V,  where  the  propeller  thrust  is  equal  to  the 
calm  water  resistance,  we  may  choose  our  coordinate  system  for 
xq,  the  surging  motion,  where  the  origin  is  also  moving  at 
velocity  V.  The  equation  for  surging  motion  can  then  be 
written  as: 

00 

(,*n)  VNio  *  "  )0  coa[~  xo+  V'+  0<“>]  ^'fa(|)r(<u)  d.u 

(6) 


The  added  mass  may  be  small  and  will  be  neglected 
lrpose  of  this  pape] 
insignificant  in  most  cases. 


for  the  purpose  of  this  paper.  The  damping  force  Nxq  will  be 


Equation  6  is  nonlinear  due  to  the  term  xQ  on  the 
right-hand  side.  This  nonlinearity  is  insignificant  in  the 
cases  of  a  ship  at  zero  speed  or  running  in  head  Beas,  but 
may  be  significant  for  a  ship  running  in  following  seas. 


Summarizing,  for  the  ship  running  in  head  Beas  or 
ip  at  zero  speed,  the 
neglected.  Then  eq.  6  becomes: 


for  the  ship  at  zero  speed,  the  terms  xQ  and  can  be 


•  M  ♦ 

C  +“  X  = 

o  m  o 


00 

L  COS  [u>(l+-^)t+  o(co)J  .  /J 0)4f2(^)  r(  cd)  dtD  (7) 

6  * 


and  the  spectrum  for  the  motion,  that  is  for  x  .  x  or  x  , 
c  000 

can  be  determined.  For  instance,  xQ  will  be: 
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(8) 


^  cos|u>(l  +  -jjr)  t  +  o 


rM 

(1+  ^)2 


doo 


If  the  same  assumptions  were  used  for  the  ship  in 
a  following  sea,  eq.  7  would  be  replaced  by: 


Xq+^o  =  j*  cos  (l-  —)  t  +  a  (oj)J  J  f 2  ij)  r((u)  dm  (9) 

0  ' 

and  instead  of  eq.  8,  the  following  surging  velocity  would 
be  found: 


r(o>) 

-  do 

(1  -  ^)2 


(10) 


Because  (l  -  equals  0  for  a  certain  circular 

o 

frequency  0  and  because will  be  a  small  value,  eq.  10 
would  produce  a  large  velocity  xq.  Moreover,  the  displace¬ 
ment  xq  computed  from  equation  9  or  equation  10  would  become 
infinite.  Therefore,  it  is  not  permissible  to  neglect  xQ 
on  the  right-hand  side  of  equation  6  for  a  ship  in  a  follow¬ 
ing  sea. 


Reference  2  shows  that  the  nonlinearity  represented 
by  the  term  xq  on  the  right-hand  side  of  eq.  6  is  significant 

in  the  case  of  the  ship  in  a  regular  following  sea  and  is 
responsible  for  the  striking  behavior  mentioned  in  the 
Introduction.  Unfortunately,  a  scientific  treatment  of  this 
behavior  in  irregular  following  seas  is  not  possible  by  the 
method  of  St.  Denis  and  Pierson.1  Therefore,  another  method 
is  proposed  which  does  not  result  in  a  complete  solution  but 
does  provide  answers  to  certain  basic  questions  about  surging 
motions  in  irregular  following  seas. 
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THE  IMPULSE  SPECTRUM 


The  way  to  such  answers  i3  opened  by  introducing 
an  impulse  spectrum  which  will  be  described  below. 

Equation  6  may  be  sufficiently  correct  to  satisfy 
the  aims  of  this  paper  provided  that  a  means  of  solving  this 
equation  can  be  found.  The  equation  is  nonlinear,  however, 
and  it  is  not  possible  to  obtain  a  solution  in  the  conven¬ 
tional  manner. 


A  different  approach  will  now  be  tried.  Let  time 
be  stopped  at  an  arbitrary  point,  t,  and  let  the  succeeding 
time  be  represented  by  t  .  Neglecting  M  ,  eq.  6  is  now 
written  as 


^o+m  Xo  =  J  cos[y  xo+cue^t+T)  +  o  U)j 7u>4  f  2(^)r(o>)  do>  ( ll) 
0 

Equation  11  can  be  written  for  a  ship  in  a  head  sea 
as  follows: 


oo 

••  n  •  (  r 

x  +—  x„  =  I  cos 
o  m  o  J  L 


^  xo+u>(l+^)  (t+r)  +o(cu) 


^o)4  f2(^;Mu))  dec 

(12) 


and  for  a  ship  in  a  following  sea: 


100 

fo>2 

COS  —  X^ 

L  g  O 


(l-|r)(t  +  T)  +0  (CD) 


)]y&7 


2(_£i)  r(cc)  dco 
(13) 


To  illustrate,  assume  that  in  a  model  test,  the 
model  initially  is  attached  to  a  carriage  moving  at  a  constant 
speed,  V,  at  which  the  mean  propeller  thrust  equals  the  mean 
resistance  (see  remarks  preceding  eq.  6).  Then,  at  an 
arbitrary  time  t,  and  r  =  0,  the  model  is  released  from  its 
attachment  to  the  carriage.  For  the  free  model,  time  is  now 
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represented  by  r  and  t  is  considered  as  a  constant.  Obvious¬ 
ly  it  is  not  possible  to  make  the  same  assumption  for  a  full 
scale  ship.  But  this  will  not  exclude  the  use  of  the  same 
method  for  the  ship,  because  in  an  irregular  sea  a  severe 
sea  can  follow  a  long  Interval  of  relatively  moderate  seas. 

Using  this  concept,  xQ,  xQ,  and  xQ  in  eqs.  11,  12, 
or  13  will  represent  functions  of  r  with  certain  conditions 
at  t  =  0.  Obviously  it  is  not  possible  to  solve  eqB.  11,  12 
and  13  in  this  form  for  the  reason  given  in  connection  with 
eq.  6.  Next,  both  sides  of  these  equations  are  Integrated 
over  r  from  t  =  0  to  r  =  T.  Equation  13  then  becomes 
T  T  oo 

J[xo+jj-  xQj  dr  Jc osj^-  xq-  <d  (l  --jp)  (t+r)+a((jD^a>4  f 2(-^)  rU )  dco 

0  0  0  ( 14) 

The  left-hand  side  of  eq.  14  becomes 

xQ  (t+T)  -  xQ(t)  +^[xo(t-n)  -  xQ(  t)]  (15) 

In  general,  the  right-hand  side  of  eq.  14  cannot 
be  integrated  and  a  complete  solution  is  not  possible.  But 
the  following  procedure  will  now  be  tried.  If  an  arbitrary 
surge  motion,  xQ,  is  assumed,  the  right-hand  side  of  eq.  14 

can  be  Integrated  with  respect  to  t.  The  result  of  this 
integration  can  then  be  compared  with  exp.  15  to  see  if  the 
assumed  xQ  will  satisfy  that  equation.  This  is  the  most  im¬ 
portant  part  of  the  proposed  method  and  the  spectrum  obtained 
as  a  result  of  integrating  the  right-hand  side  of  eq.14  with 
respect  to  r  is  called  an  impulse  spectrum. 

The  most  interesting  question  about  surging  motion 
in  a  following  sea  is  whether  the  ship  can  be  accelerated  by 
the  waves  and  maintain  a  speed  higher  than  its  mean  speed  for 
an  extended  time.  To  examine  this,  first  assume  an 


dr 
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accelerated  motion,  xQ,  for  computing  the  Impulse  spectrum 
and  then  test  this  assumption  to  see  if  such  a  motion  satis¬ 
fies  eq.  14. 

Now  for  x  ,  assume  an  accelerated  motion  with  a 
constant  acceleration,  b. 


x  =  b  r  2 
o  2 


(16) 


Then,  the  right-hand  side  of  eq.  14  becomes: 


Re 


oo 

I 

0 


i  [  -0)  (l~  )  t+  0  (ai 


y 


cu4f  2  (-j-)  r(o>  )dco 


(i 

Jo 


Now  the  Integration  over  r  can  be  performed: 
T 


/ 


i[-(i-f)r+f  i  t-; 


dT 


If  b  =  0,  expression  18  becomes: 

,oT  / ,  oVi 


•if  fi-f) 


0).  ( ,  cbV\ 
2 


t( ,  ciV\  .a)2  b 

■“f1— >T+r?T 

(17) 

(18) 

( 19a) 


i 


If  b/o,  expression  18  becomes: 


01 


eiudu 


(19b) 


■d-^v) 


/jsT 


g  '  V2b 


•  The  Integral  in  exp,  19b  is  known  as  a  Fresnel 

Integral  and  is  tabulated  in  ref.  4.  If  the  computations  are 

to  be  performed  on  an  electronic  digital  computer,  it  will  be 

iu2 

more  convenient  to  integrate  after  e  is  expanded  in  a  series. 
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The  exponential  function  in  exps.  19a  and  19b  can  be 


included  in  a  new  phase,  a(o>),  of  the  resulting  representa¬ 
tion  of  the  impulse.  Without  this  exponential  function, 
exp. 19b  is  T  times  a  dimensionless  complex  function  of 


Then  exp.  17  can  be  replaced  by 


The  phase  0  in  exp.  21  is  not  the  same  as  that 
in  exp. 17  but  this  is  not  important  since  in  each  case  a  1b 
a  random  value.  The  spectrum 


co4f2(-£)  T2  (iz+y2)  r(a>)  (22) 

may  be  called  the  impulse  spectrum  for  the  surging  motion. 

The  impulse  (exp.  2l)is  represented  as  an  arbitrary 
superposition  of  a  finite  or  infinite  number  of  harmonic  com¬ 
ponents  in  the  same  manner  as  the  irregular  seaway  itself. 

It  is  dependent  on 

1.  the  seaway  spectrum,  r 

2.  the  forcing  function,  f  (-^) 

3.  the  time  for  which  the  impulse  is  computed,  T 

4.  the  assumed  surge  motion,  xq 
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Now  eq.  14  must  be  satisfied.  Using  the  assumed 
motion  of  eq.  16  on  the  left-hand  side  of  eq.  14  as  repre¬ 
sented  in  exp.  15: 


bT+ 


f°0 

cos 

0 


f2(_)T2(  I2+y2  ) 


r(o))  dew 

(23) 


This  equation  is  satisfied  only  for  a  unique  value 
of  b  and  it  can  be  used  to  estimate  this  acceleration,  b,  as 
a  function  of  the  time  Interval,  T.  The  influence  of  the 
nonlinearity  due  to  the  term  xq  on  the  right-hand  side  of 
eq.  14  will  be  Included  in  this  estimate. 


APPLICATION  TO  THE  BEHAVIOR  OP  A  SHIP 
IN  A  REGULAR  FOLLOWING  SEA _ 

To  demonstrate  what  can  be  expected  of  the  proposed 
method,  first  the  case  in  regular  following  waves  is  treated. 
For  this  case,  a  result  found  by  a  more  exact  method  is 
known.2  For  this  special  case,  eq.  23  becomes: 

bT  +  mlT2=  cos  [a>(l-|¥)t  +  o(.)]  u>2f  (-|)  t/i^+Y2  h 

because  in  this  case  do>  must  be  replaced  by  the  wave 

amplitude  h. 

0>2f 

value  b  represents  the  amplitude  of  the  acceleration  which 
can  be  produced  by  the  wave.  This  can  be  seen  from  eqs .  2 
and  5.  Then  eq.  24  is  rewritten: 

bT+  ^  |  Ta=  b  T  ^/^HhY2  cos  [  o>(  1~  )  t  +  o(u>)J  (25) 

Remember  the  meaning  of  b:  it  represents  an 
average  acceleration  over  a  time  interval  T  which  begins  at 
the  time  t. 


by  b.  This 


It  is  convenient  to  replace 


[ 
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We  are  most  interested  in  the  maximum  value  of  the 
impulse.  For  this  maximum  value,  the  harmonic  function  in 
eq.  25  is  replaced  by  1.  This  harmonic  function  only  serves 
to  show  that  the  impulse  depends  on  the  time  t  at  which  the 
model  is  released  from  the  carriage  and  that  this  dependence 
is  harmonic.  For  this  most  interesting  case,  eq.  25  becomes: 

bT+—  |  T2  =  bT  aJiz+Y3  ( 26) 

Remember  that  I  and  Y  are  functions  of  a IT,  (l  -~) , 

C> 

—  .  Now  for  the  regular  waves,  given  by  and  5,  and  fob 

o 

a  mean  velocity  V,  a  time  interval  T  is  chosen  and  the  maximum 
value  of  the  average  acceleration  b  for  this  time  interval 
or  the  Increment  of  velocity  bT  is  determined  from  eq.  26. 
This  computation  is  carried  out  for  different  time  intervals 
T  and  the  results  are  presented  in  Fig.  2. 

In  Fig.  2  the  increment  of  velocity  produced  by  the 
waves  in  the  time  interval  T  1b  plotted  versus  T.  This 
diagram  shows  that  for  small  values  of  T,  the  increment  of 
velocity  does  not  depend  on  the  velocity  V.  In  this  range 
the  average  acceleration  b  equals  the  maximum  possible 
acceleration  b.  But  for  larger  time  intervals  T  the  average 
acceleration  b  is  smaller  than  b  because  the  location  of  the 
model  relative  to  the  wave  is  changed  within  this  interval  T. 
The  influence  of  the  velocity  V  on  the  increment  of  velocity 
is  very  significant.  If  the  velocity  V  equals  the  wave 

velocity  C=§,  the  parameter  on  the  curves  equals  0  and  the 

Increment  of  velocity  is  relatively  small.  But  the  higher 
velocity  will  continue  for  a  relatively  long  time.  If  the 
velocity  V  is  smaller  than  the  wave  velocity,  the  parameter 
on  the  curves  increases  and  the  increment  of  velocity  also 
increases.  But  if  the  velocity  V  becomes  too  small,  a  very 
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striking  change  can  be  seen.  Then  the  increment  of  velocity 
remains  small,  and  furthermore,  this  increment  exists  for 
only  a  very  short  time.  It  can  be  concluded  from  Fig.  2 


that  for  values  of  the  parameter  j^(  smaller  than 

1.1,  the  model  will  be  accelerated  to  such  an  extent  that  it 
will  attain  the  wave  velocity  and  run  with  the  wave.  But 

for  parameters 


larger  than  1.1,  the  increment 


of  velocity  remains  small  and  the  model  can  never  maintain 
the  wave  velocity. 


Nearly  the  same  result  is  obtained  in  ref.  2. 
Quantitatively  the  agreement  is  not  perfect.  But  a  perfect 
agreement  could  not  be  expected  because  the  question  is  not 
exactly  the  same  as  in  ref.  2.  Reference  2  poses  the  ques¬ 
tion  as  to  the  condition  under  which  a  wave  will  force  a  ship 
to  run  at  the  wave  velocity.  On  the  other  hand,  the  method 
developed  here  is  used  to  find  what  maximum  increment  of 
velocity  can  be  imposed  on  a  ship  by  a  wave  in  a  given  time 
interval.  In  spite  of  this  difference,  comparison  of  the 
results  shows  that  the  method  developed  in  this  report  is 
reasonable  and  can  be  used  to  obtain  useful  answers  for  a 
model  or  a  ship  in  regular  following  waves. 

Of  course,  the  increment  of  velocity  may  be  very 
small  and  may  even  be  negative,  depending  upon  the  time  t  at 
which  the  model  is  released.  But  we  are  primarily  interested 
in  the  largest  possible  increment  of  velocity  and  therefore 
have  replaced  the  harmonic  function  in  eq.  25  with  a  value  of 
unity. 

Strictly  speaking,  due  to  the  assumption  of  eq.  16, 
Fig.  2  should  be  used  only  in  the  region  where  the  accelera¬ 
tion  is  constant  (close  to  the  line  b  =  b) .  However,  a  good 
approximation  is  obtained  if  each  curve's  use  is  restricted 
to  the  range  below  the  point  of  maximum  increment  of  velocity. 
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It  should  be  remembered  that  nonlinearities  have 
an  important  effect  on  these  results.  An  accelerated  ship 
remains  on  the  slope  of  a  wave  which  produces  the  accelera¬ 
tion  a  longer  interval  of  time  than  does  a  non-accelerated 
ship.  This  effect  is  included  in  the  computation  only  if  the 
nonlinearities  are  taken  into  account. 

APPLICATION  OP  THE  PROPOSED  METHOD  TO  CASES  IN  WHICH 
A  LINEARIZATION  MAY  BE  ALLOWED 


For  the  ship  running  in  head  seas  or  for  the  ship 
at  zero  speed  the  nonlinearity  is  not  significant  and  there¬ 
fore  can  be  neglected.  Now  it  will  be  demonstrated  that  the 
proposed  method  delivers  useful  results  for  this  case. 


For  the  model  in  regular  waves,  eq.  25  or  eq.  26 
must  be  used.  But  for  cases  in  which  the  nonlinearity  can  be 
neglected,  the  more  simple  expression  for  (I  +  iY)  written 
in  exp. 19a  can  be  used  Instead  of  the  more  complicated  ex¬ 
pression  written  inexp.  19b: 


(I  +  1Y) 


sln[f(l  +  £?)] 

jcdT  ,  <jjV\  for  head  seas  (27) 

2  v±  +  g  ' 


Then  eq.  26  for  determining  the  maximum  value  of  the  average 
acceleration  b  or  of  the  increment  of  velocity  (bT)  becomes: 


bT  +  —  T2  =  bT 
m  d 


(28) 


N 

The  Influence  of  the  damping  coefficient  —  is  small  if 
application  is  restricted  to  the  range  below  the  first  maxi¬ 
mum.  The  maximum  possible  increment  (for  ~  =  0)  amounts  to: 
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(29a) 


for  (bT) 


max 


(29b) 


Solution  of  the  linearized  equation  of  surging  motion  with 
zero  damping  force  in  a  regular  head  sea  gives: 

1.  an  amplitude  of  oscillatory  surging  velocity 
one-half  the  value  obtained  from  eq.  29a,  and 

2.  a  period  double  the  value  obtained  from 

eq.  29b. 

This  confirms  the  relations  for  a  harmonic  oscil¬ 
lation,  where  the  maximum  increment  of  velocity  equals  two 
times  the  velocity  amplitude,  and  it  occurs  during  one-half 
of  the  period  of  oscillation,  as  shown  in  Pig.  3» 

For  the  case  in  an  irregular  head  sea,  the  spectrum 
on  the  right-hand  side  of  eq.  23  becomes,  without  the  non¬ 
linearity. 


This  result,  for  a  vanishing  damping  force,  is  in 
complete  agreement  with  the  result  obtained  in  the  normal  way 
from  the  linearized  equation  of  motion  (eq.  7)-  Also,  for 
the  case  of  a  vanishing  damping  force,  exp.  30  represents  the 
spectrum  of  the  increment  of  velocity.  On  the  other  hand, 
eq.  8  is  the  result  for  the  velocity  of  a  straightforward 
solution  of  eq.  7.  The  increment  in  any  time  interval  T, 
therefore,  can  be  determined  for  this  linearized  case  direct¬ 
ly  from  eq.  8. 


R-929 

-13- 


Increment  of  velocity  =  xQ  ( t  +  T)  -  x  (t)  = 


=  Jjcos 
0  ' 


[o>(  l+|Y)  ( t+T)+a(u>)]  -cos[cj(  l+sp)  t+  a 

■J  (1  +f)* 


v’O) 


(31) 


do 


The  spectrum  of  the  increment  of  velocity  found  in 
this  normal  way  is  the  same  as  the  impulse  spectrum  (  exp.  30) . 

For  the  cases  discussed  in  this  section  and  the 
previous  one,  there  is  no  advantage  in  using  the  proposed 
method  in  place  of  the  method  developed  in  refs.  1  or  2. 
However,  the  former  method  can  not  be  used  in  treating  the 
surging  motion  in  an  irregular  following  sea  since  the  in¬ 
fluence  of  nonlinearity  is  not  included. 


IMPULSE  FOR  ANOTHER  ASSUMPTION  OF  THE  VELOCITY  x 
- — — — — - o 

In  the  section  entitled  "impulse  Spectrum,"  a 
motion  with  a  constant  acceleration  b  within  the  time  in¬ 
terval  T  is  assumed  and  then  the  impulse  is  calculated.  It 
is  known  from  refs,  2  and  3  that  In  a  regular  following  sea 
the  ship  is  accelerated  only  to  the  wave  speed  and  after  that, 
it  runs  with  the  wave  for  an  unlimited  time.  Figure  2  also 
shows  a  similar  trend. 

Therefore,  it  is  also  Interesting  to  calculate  the 
impulse  for  a  velocity  xQ  which  remains  constant  after  a 
certain  increment  of  velocity  is  reached. 
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For  this  case,  exp .  l8  must  be  supplemented  by  a  second  term: 


fTl  i[-oj(l-|?)T+  ^  §r2]  *2  i[-a)(l-|Z)T+  ^iT^bT^r-Tj] 

I  e  dT+  \  e  dT 

*  ~  / _ \ 


Expression  33  then  takes  the  place  of  exp.  19b  to  become 
the  following: 


It  may  be  supposed  that  the  assumption  of  an  interval  of 
constant  velocity  following  acceleration  will  lead  to  a  pre¬ 
diction  of  a  ship  running  at  the  dominant  wave  speed  for  an 
extended  period  of  time.  Instead  of  exp.  22,  the  resulting 
Impulse  spectrum  becomes: 

cu4  f2&)  T  2  (loa  +  Y  a)  rU)  (35) 

u  i  2  a 
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where  T^2  1 22  +  Yg2  J  represents  the  square  of  the  absolute 

value  of  exp.  3^. 

Discussion  of  the  irregular  sea  case  will  be  con¬ 
tinued  later;  the  following  applies  only  to  a  regular  sea. 

It  is  assumed  that  the  ship  is  accelerated  to  the  wave  speed 
and  the  question  is  whether  or  not  the  ship  will  continue  to 

run  at  the  wave  speed,  (wave  speed  =  ■§) 

bT.  =  i.  -  V  (36) 

1  OD 

This  equation,  put  in  exp.  3^,  results  in: 


(37) 


This  means  that  during  the  second  time  Interval  Ti< 
t  <  T2  ,  the  Impulse  grows  linearly  with  the  time,  or  in  other 
words,  the  force  in  the  longitudinal  direction  remains  con¬ 
stant.  Therefore,  if  the  condition  of  eq.  26  can  be  ful¬ 
filled  for  the  first  time  interval,  the  ship  can  run  with 
the  wave  an  infinite  time.  Also  this  result  is  very  close 
to  the  results  obtained  in  refs.  2  and  3* 

A  mean  wave  speed  will  be  dominant  in  an  irregular 
sea.  If  the  model  can  be  accelerated  by  the  sea  to  this  mean 
wave  speed,  a  long  run  at  this  wave  speed  appears  possible. 
Therefore,  it  can  be  expected  that  the  approach  developed  in 
this  section  will  be  useful  for  the  case  in  an  irregular  sea. 


SOME  RESULTS  FOR  NONLINEAR  BEHAVIOR  IN 
AN  IRREGULAR  FOLLOWING  SEA _ 


In  the  fourth  and  fifth  sections  it  was  shown  that 
the  proposed  method  leads  to  useful  results  for  the  oscilla¬ 
tory  surging  motion  in  an  irregular  head  sea  as  well  as  for 
the  striking  nonlinear  behavior  of  a  model  in  regular  follow¬ 
ing  waves.  Now  the  method  will  be  used  to  investigate  the 
nonlinear  behavior  in  an  irregular  following  sea. 

The  numerical  computations  are  carried  out  on  the 
IBM  1620  electronic  computer  at  Stevens  Institute  of 
Technology. 

A  two-dimensional  modified  Neumann  spectrum  is 
assumed  for  the  irregular  seal 

2gj~ 

rM  e'“<“  (38) 

The  constant  C  amounts  to  32  ft2  sec-5  for  the  Neumann 
spectrum,  but  for  the  computations  performed  for  this  study, 

C  was  reduced  by  the  factor  0.75  to  24  ft2  sec-5.  The  in¬ 
fluence  of  the  three-dimensional  nature  of  the  real  sea  on 
the  surging  force  is  estimated  by  this  reduction. 

Equation  3  is  used  for  the  surging  force  function. 
Then  the  impulse  spectrum  amounts  to  ( see  exp.  22) : 


where  T2(l2  +  Y2)  represents  the  square  of  the  absolute  value 
of  the  Integral  ( exp.  18) .  It  can  be  taken  as  in: 
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1.  exp.  19a,  if  nonlinearity  is  excluded, 

2.  exp.  19b,  if  nonlinearity  is  included  and 
eq.  16  is  used  for  the  assumed  motion, 

3.  exp.  34,  if  the  nonlinearity  is  included  and 
eq.  32  is  used  for  the  assumed  motion. 

The  area  under  the  impulse  spectrum  is  determined  by: 


00 


This  area  depends  on  the  wind  speed  u,  ship  speed  V,  ship 
length  L  and  the  assumed  motion.  The  root  of  this  area  is 
written  as: 


where  C  =  the  reduced  factor  for  the  Neumann- spectrum  = 

24  ft2  sec"5 

L  =  the  length  of  the  ship  in  feet 
g  =  the  gravitational  acceleration 
D  =  a  dimensionless  number 

D  is  computed  from  the  Impulse  spectrum  as  a  function  of  the 
assumed  motion, of  the  Froude  number  v/  ^  gL  for  the  speed 
of  the  ship  and  of  the  dimensionless  parameter  u/  ^  gL  for 
the  wind  speed. 

For  the  assumed  motion,  the  following  must  be 

chosen: 

In  Case  1.  only  the  time  interval  T  (see  exp. 19a) .  (This 

case  is  computed  only  as  a  basis  for  comparison.) 

In  Case  2.  the  time  interval  T  and  the  acceleration  b 
(see  exp. 19b) . 
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In  Case  3.  the  time  interval  T1#  the  acceleration  b  and  the 
time  Interval  T2  (see  exp.  34). 

The  impulse  is  a  random  value.  If  a  Rayleigh  dis¬ 
tribution  is  assumed,  the  amplitude  of  the  impulse  therefore 
amounts  to: 

C1/2(L)3/4D  (42) 

g 

where  a  =  0.886  for  the  mean  value  of  the  amplitudes 

a  =  1 . 4l  for  the  average  of  the  l/3  highest  amplitudes 

a  =  1.80  for  the  average  of  the  l/lO  highest  amplitudes 

a  =  2.25  for  the  average  of  the  l/lJO  highest 

ampli  tudes 

In  each  case,  the  impulse  must  equal  the  left-hand 
side  of  eq.  23.  The  damping  term  in  this  equation  may  be 
neglected  for  many  cases.  In  such  cases,  it  follows  from 


eq.  23  and  exp.  42  that: 


/  T  3/ 4 

bT  =  aCl/2(|)  D 

(43) 

V*. 

or 

aCl/aL1/4 

(44) 

If  the  nonlinearity  Is  excluded  as  in  Case  1, 

the 

value  D  is  independent  of  the  acceleration  b.  Then  this 
value  D  can  be  computed  as  a  function  of  the  Froude  number, 
the  dimensionless  wind  speed  and  the  time  interval  T. 

Figure  4  shows  the  results  of  such  a  computation.  The  in¬ 
crement  of  velocity  (bT)  can  be  determined  from  this  figure 
for  any  ship's  length  L  and  for  any  value,  "a."  If  the  non¬ 
linearity  is  included  as  in  Case  2,  the  left-hand  side  of 
eq.  44  must  be  determined  first.  Only  the  left-hand  side  of 
this  equation  depends  on  the  assumed  motion.  Fig.  5  shows  the 
result  of  such  a  computation  for  Case  2.  Fig. 5  Is  valid  only 
for  Fr**0 . 2/ and  u2/gL=l/3.  But  the  assumed  motion  is 
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varied  by  varying  the  parameters  and  T  tj b/L.  Now  it 

can  be  determined  from  this  figure  which  assumed  motion  is 
possible  for  a  given  right-hand  side  of  eq.  44.  For  instance, 
if  a  ship  of  L  =  220  ft  is  given  and  2.25  is  chosen  as  factor 
"a,"  the  value  of  the  right-hand  side  of  eq.  44  amounts  to: 


_ 325/4 _ 

2.25  x  241/2  x  2201/4 


1.79 


All  points  in  Fig.  5  having  an  ordinate  of  1.79  represent 
for  the  given  ship,  Froude  number  and  wind  speed,  an  average 
of  the  l/lOO -  highest  motions.  If  1.80  is  chosen  as  the 
value  of  the  factor  "a,"  the  right-hand  side  of  eq.  44  amounts 
to  2.24  and  all  points  in  Fig.  5  having  such  an  ordinate 
represent  an  average  of  the  l/lO-highest  motions.  Figure  6 
shows  both  of  these  motions  and  it  also  shows  the  same 
motions  determined  without  the  influence  of  the  nonlinearity 
and  computed  from  Fig.  4.  The  same  procedure  is  carried  out 
for  other  wind  speeds  and  Froude  numbers.  Figure  7  shows 
results  analogous  to  Fig.  6  but  for  a  more  severe  sea. 

Figures  6  and  7  illustrate  that  the  influence  of 
nonlinearity  is  stronger,  the  more  severe  the  sea  is  or  the 
larger  the  factor  "a"  is.  The  average  of  the  l/lOO-highest 
motions  will  be  a  larger  motion  than  the  average  of  the 
l/lO -highest  motions  and  the  Influence  of  the  nonlinearity 
must  be  stronger  for  the  larger  motion.  Therefore,  the 
ratio  of  the  average  of  the  l/lOO-highest  motions  to  the 
average  of  the  l/10-hlghest  motions  is  larger  than  the  ratio  of 
the  "a"  values  (2.25/1.80  =  1.25)  for  the  nonlinear  motion. 
Figure  8  shows  results  for  other  wind  speeds. 

The  damping  term  on  the  left-hand  side  of  eq.  23  is 
neglected  in  computations  for  Figs.  4  through  8.  The  motion 
will  be  reduced  by  this  term.  It  is  supposed  that  this 
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influence  is  unimportant  for  u2/gL  =  l/5  and  l/4  in  Pig.  8, 

but  may  be  Important  for  u2/gL  =  l/3  and  l/2  in  this  figure. 

However,  this  influence  will  be  discussed  later.  Now  it  may 
be  concluded  only  that  for  a  Proude  number  of  0.2/ ^[2~,  the 
increment  of  the  velocity  becomes  very  large  for  wind  speeds 
of  u2/gL  =  l/3  or  more. 

In  the  next  step  Case  3  is  treated,  that  is,  eq.  32 
is  taken  for  the  assumed  motion.  First  the  damping  term  on 
the  left-hand  side  of  eq.  23  is  neglected.  Then  eqs.  43  and 
44  must  be  fulfilled  for  this  case.  The  difference  compared 

to  Case  2  is  only  that  T  in  these  equations  must  be  replaced 

by  T-^  and  that  the  assumed  motion  is  expressed  in  D  by  the 
parameters  b,  Tx  and  T2 .  It  appears  best  to  use  the  results 
already  obtained  for  Case  2. 

In  Case  2  a  probable  motion  xq  =  br2/2  or  xq  =  br 

was  obtained  in  the  time  interval  T  for  given  values  of 
the  Froude  number  Fr,  the  wind  speed  u,  the  length  of  the 
ship  L,  and  the  value  "a." 

It  is  now  asked  what  the  impulse  is,  if  it  is 
assumed  that  for  the  same  pre-conditions,  after  the  time 
interval  r  =  T^,  the  ship  runs  with  a  constant  velocity 
(V  +  bT^  .  The  following  form  of  eqs.  43  or  44  is  used  to 
represent  the  results  of  such  a  computation: 


,  T  l/4 

aCl/2  L_  D 


S 


5/4 


(45) 


The  right-hand  side  of  this  equation  is  computed 
from  the  impulse  spectrum  and  it  is  plotted  versus  the  time 
T2  VSA  m  Pig.  9*  The  left-hand  side  of  this  equation  is 
determined  from  the  assumed  motion  and  plotted  as  dashed 
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lines  in  this  figure.  Equation  45  is  fulfilled  only  if  the 
corresponding  curves  in  the  figure  agree. 

This  agreement  may  be  close  enough  for  the  first 
time  interval.  The  results  for  this  time  interval  do  not 
agree  exactly  with  the  corresponding  results  represented  in 
Pigs.  6  to  8.  The  reason  for  this  discrepancy  is  only  that 
to  save  time  the  computations  are  not  repeated  with  the 
parameters  b  and  T  which  can  be  read  from  Pig.  8. 

However,  the  agreement  of  corresponding  curves  in 
Pig.  9  is  not  close  enough  for  the  time  r  >  T-^. 

In  spite  of  these  discrepancies,  it  can  be  con¬ 
cluded  that  for  the  Proude  number  Fr  =  0.2/  ^/~2",  the  ship's 

length  L  =  220  ft  and  for  wind  speeds  of  ua/gL  =  l/5  and  l/4, 
the  ship  will  run  with  a  higher  speed  only  a  short  time.  On 
the  other  hand,  it  is  possible  that  for  u2/gL  =  l/3  and  l/2, 
the  real  motion  lies  somewhere  between  the  computed  and  the 
assumed  motion,  and  also  that  the  ship  will  run  with  a  higher 
speed  a  long  time. 

It  may  be  mentioned  at  this  point  that  the  curves 
for  the  computed  motion  and  for  u2/gL  =  l/3  and  l/2  in  Pig.  9 
for  a  time  r  steadily  Increasing  will  reach  a  maximum  and 
after  that  will  descend.  However,  further  discussion  should 
include  the  Influence  of  the  damping  term  in  eq.  23. 

Furthermore,  it  may  be  mentioned  that  the  assump¬ 
tion  for  Case  3  seems  to  lead  to  very  significant  results 
for  behavior  in  an  irregular  sea. 
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INFLUENCE  OF  SHIP'S  RESISTANCE  ON  NONLINEAR  BEHAVIOR 


The  damping  term  on  the  left-hand  side  of  eq.  23 
is  due  to  the  resistance  of  the  ship.  This  term  is  neglected 
in  the  preceding  numerical  computation.  But  such  an  omission 
is  not  allowable  for  cases  in  which  the  ship  runs  at  a  high 
speed  for  a  long  time. 

Moreover,  the  increment  of  velocity  may  be  very 
large  (Fig.  9)  and  therefore  the  linearized  form  of  the 
damping  term  used  in  eq.  23  may  not  be  correct  for  such 
large  increments  of  velocity. 

For  estimating  the  influence  of  damping,  a  quad¬ 
ratic  law  is  assumed  for  the  resistance.  The  damping  term 
on  the  left-hand  side  of  eq.  23  represents  the  difference 
between  resistance  R  and  propeller  thrust  T. 

For  the  following  computation  it  is  assumed 

R-T  =  —  f  (V  +  xj2  -  V2]  (46) 

y2  l-  O  -I 


where  Rv  =  a  constant  value.  This  value  will  not  equal  the 
resistance  of  the  ship  for  the  speed  V  because  the  ship's 
propeller  will  run  with  about  a  constant  number  of  revolu¬ 
tions  even  when  the  ship  is  accelerated  by  the  waves.  The 
value  Rv  in  eq.  46  will  therefore  be  larger  than  the  resist¬ 
ance  for  the  speed  V,  and  it  must  be  determined  by  a  resist¬ 
ance  test  of  a  model  equipped  with  a  propeller  running  with 
about  a  constant  number  of  revolutions. 


The  left-hand  side  of  eq.  13  is  replaced  by: 


•  • 


(2Vxr 


(47) 
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For  Case  3>  which  has  already  been  treated,  a 
motion  determined  by  eq.  32  is  assumed: 


Xo  =  br 

for 

0  SrsT 

i0  -  Mj 

for 

T^rsT, 

The  integration  of  eq.  47,  similar  to  the  integra¬ 
tion  of  the  left-hand  side  of  eq.  14,  leads  to: 


1 X  +— ^— ( 2Vx+x2)"] 
J  L  0  mV2  °  0  J 


dr  =  bT, 


1+  -L-  |2V(T2-it)  +bT1(T2-|ri)] 

(48) 


The  left-hand  side  of  eq.  23  will  therefore  be  re¬ 
placed  by  eq.  48.  The  frequency  co  does  not  appear  in  this 
expression  and  therefore  it  is  only  necessary  in  eq.  43  to 
45  to  divide  D  by 


or,  represented  by  nondlmensional  values: 


(49) 


(50) 

The  new  parameter  — — —  is  introduced.  Unfortu- 

mgFr2 

nately  the  magnitude  of  this  parameter  will  vary  for  different 
ships.  Therefore,  the  following  numerical  results  will  show 
only  a  general  trend,  and  these  results  should  not  be  used 
for  any  actual  ship.  Rv/mgFr2=  0.15  is  used  for  the  follow¬ 
ing  numerical  computations .  This  may  be  a  large  value,  but 
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it  must  be  remembered  that  the  influence  of  the  propeller  is 
included  in  this  value. 

Figure  10  is  obtained  by  dividing  the  computed  val¬ 
ues  plotted  in  Fig. 9  by  the  factor  from  exp.  50.  A  comparison  of 
the  two  figures  shows  that  the  influence  of  the  damping  term 
is  relatively  small  during  the  first  time  interval  but  that 
it  is  of  great  importance  for  the  next  time  interval.  Again 
it  can  be  seen  from  Fig.  10  that  in  a  mild  sea,  a  high  speed 
is  possible  only  for  a  short  time.  It  should  be  considered 
accidental  that  for  ua/gL  =  l/3  the  computed  motion  is  so 
close  to  the  assumed  motion.  If  the  increment  of  velocity 
for  u2/gL  =  1/3  and  1/2  in  Fig.  10  appears  unrealistically 
high,  it  should  be  remembered  that  analogous  results  were 
presented  in  refs .  2  and  3  for  the  case  in  regular  waves . 

Also  in  ref.  2,  a  very  high  increment  of  velocity  was  found 
for  a  severe  sea  and  even  in  ref.  3  a  relatively  high  incre¬ 
ment  was  found,  even  though  a  mild  sea  was  used  in  the  ex¬ 
periments.  For  instance,  in  Fig.  14  of  ref.  3  an  increment 
of  bT/^r  =0.l66  is  found  for  a  wave  speed  of  8.9  ft/sec 
although  the  ratio  of  wave  height  to  wave  length  was  only 
1:  41. 


AN  ESTIMATE  OF  THE  MAXIMUM  DURATION  OF  A  RUN 
WITH  A  HIGHER  SPEED  THAN  THE  CALM  WATER  SPEED 

Figure  10  indicates  that  when  u2/gL  =  1/3  and  l/2, 
it  is  possible  for  the  ship  to  maintain  a  high  speed  for  a 
long  Interval.  For  cases  in  which  the  second  time  interval 
(T2  -  Tj.)  is  very  long,  it  appears  possible  to  make  an 
estimate  of  this  time  interval  by  a  simple  method. 

Using  the  case  of  u2/gL  =  l/2  in  Fig.  10  as  an 
example,  the  problem  is  to  determine  the  time  at  which  the 
curves  of  the  computed  motion  and  the  assumed  motion  will 
intersect.  The  method  used  in  this  section  is  not  expected 
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to  result  in  a  very  accurate  estimate,  but  a  rough  approxima¬ 
tion  only.  Nevertheless,  the  results  should  help  broaden 
our  understanding  of  the  nonlinear  behavior  of  a  ship. 


The  method  of  estimation  is  based  on  the  fact  that 
the  factor  Tx2  ( I22  +  Y22)  in  the  impulse  spectrum  for  a 
long  time  (t2  -  Ti)  can  be  approximated  by  (see  exps.  34  and 
35): 


T12  (l22  +  Y22) 


sin  [|  (T2  -  Ti)  (-1  +  + 


u/Ti 


t>I 


£  ("I  +  ~  +“>T1  |) 


(51) 


This  is  because  for 


(-1  +  |Y  +CUT!  |)->  0  (52) 

the  part  of  T12  (l22  +  Y22)  represented  by  eq.  51  becomes  a 
very  large  value.  Although  this  part  of  T12  (l22  +  Y22)  de¬ 
creases  quickly,  if  a>  becomes  smaller  or  larger  as  is  needed 
for  exp.  52,  this  part  will  be  the  dominant  part  in  the  cases 
mentioned.  The  integrand  in  exp.  40  will  be  relatively  small 
for  all  values  of  os  except  a  very  small  range  near 


1 

“o  “  V  7  TI  b 
S+  g 


(53) 


For  such  a  condition  it  will  be  allowable  to  re¬ 
place  exp.  40  by: 


The  integration  is  carried  out  and  the  following 
is  obtained: 


The  factor  from  exp. 50  for  the  left-hand  side  of 
eq.  23  is  approximated  for  a  very  large  value  of  Ta  by 


Now,  both  of  these  approximations  are  used  to 
establish  the  following  equation  in  place  of  eq.  43: 


The  time  interval  (t2  -  Ti)  follows  from  that 

equation: 


(Ta-Tx)  « 


(58) 
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These  results  are  plotted  In  Pig.  11  versus (co2L/2g) . 
In  the  upper  part  of  this  figure,  the  increment  of  velocity- 
computed  from  eq.  53  is  plotted  and  in  the  lower  part  the 
time  computed  from  eq.  58  is  plotted. 


The  following  three  steps  must  be  made  in  order  to 
understand  Pig.  11; 

1.  It  must  be  determined  to  which  increment  of 
velocity  the  ship  can  be  accelerated  by  the  sea.  Such  a 
determination  will  lead  to  results  like  those  plotted  in 
Pigs.  5  to  8. 

2 .  The  value  0>2L/2g  must  be  determined  from  the 

upper  part  of  Pig.  11  for  the  increment  of  velocity  just 
found . 


3.  An  estimate  of  the  time  for  the  run  of  the  ship 
at  increased  speed  can  be  obtained  for  this  value  of  u>aL/2g 
from  the  lower  part  of  Pig .  11 .  0 


For  instance,  for  the  case  u2/gL  =  l/2  represented 
in  Pig.  10,  the  increment  of  velocity  produced  in  the  first 
interval  amounts  to  0.35.  For  this  increment  the  value  2.04 
for  o>2L/2g  is  obtained  from  the  upper  part  of  Fig.  11  and 
for  the  same  abscissa,  the  time  (Tz  -  Ti)  ^g/L  =  73-5  is 
determined.  Therefore,  both  curves  plotted  in  Pig.  10  for 
ua/gL  =  l/2  will  intersect  after  the  very  long  time  of  73.5 
and  the  real  motion  may  lie  anywhere  between  the  two  plotted 
motions . 


The  same  procedure  delivers  a  time  (Ta  -  Tx 
for  the  case  u2/gL  =  l/3  represented  in  Fig.  10. 


For  the  other  cases,  u2/gL  =  l/4  and  l/5,  this 
estimate  is  not  possible  because  the  increment  of  velocity 
produced  in  the  first  time  interval  is  too  small  and  lies 
outside  of  Pig.  11. 


Figure  11  may  also  be  used  to  find  the  increment  of 
velocity  which  will  result  in  the  longest  run  at  high  speed. 
The  probability  of  such  a  long  run  occurring  is  somewhat 
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smaller  than  that  Implied  by  the  factor  "a."  This  Is  because 
the  equation  used  In  this  computation  is  subject  only  to  the 
condition  that  the  assumed  motion  and  the  computed  motion 
coincide  at  the  end  points  of  the  time  interval.  The  proba¬ 
bility  could  be  maintained  only  if  the  assumed  and  computed 
motions  were  the  same  over  the  entire  time  interval .  On  the 
other  hand,  it  might  be  possible  to  estimate  the  decrease  in 
the  probability  of  obtaining  a  long  run,  but  it  is  believed 
that  the  results  of  such  an  estimate  would  not  be  of  great 
importance . 

CONCLUSIONS 

1.  It  appears  impossible  to  determine  the  surging 
motion  in  a  severe  following  sea  without  including  the  non¬ 
linearity  in  such  a  determination. 

2.  An  enormous  change  in  the  trend  of  the  surging  motion 
due  to  the  nonlinearity  seems  possible. 

3.  A  large  increase  in  speed  and  a  run  at  this  higher 
speed  over  a  relatively  long  Interval  of  time  can  be  pro¬ 
duced  by  a  following  sea  on  a  ship  running  slower  than  the 
dominant  waves. 

4.  The  increment  of  speed  and  the  duration  of  the  run 
at  higher  speed  can  be  estimated  by  the  method  described  here. 

5.  The  surging  motion  itself  is  not  dangerous  for  the 
ship.  But  it  is  believed  that  the  nonlinear  surging  motion 
treated  here  may  be  connected  with  dangerous  motions  in  other 
modes.  It  is  believed  that  the  treatment  of  surging  motion 
presented  here  will  provide  a  valuable  starting  point  for 
future  investigations  of  other  more  dangerous  motions. 
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BROACHING  OP  A  SHIP  IN  A  SEVERE 
IRREGULAR  FOLLOWING  OR  QUARTERING  SEA 


INTRODUCTION 

Broaching  of  ships  in  severe  following  or  quarter¬ 
ing  seas  has  been  observed  and  examples  of  broaching  are 
described  in  ref.  3-  Broaching  of  a  ship  is  very  dangerous 
and  must  be  avoided.  Therefore  it  is  important  to  know  the 
conditions  which  lead  to  such  an  occurrence. 

The  problem  was  discussed  in  connection  with  non¬ 
linear  surging  motion  in  refs.  2  and  3*  However,  broaching 
is  a  very  complicated  problem  and  not  even  an  approximate 
solution  is  known. 

Some  aspects  of  the  problem  are  discussed  in  this 
part  of  the  report.  Neither  a  complete  solution  nor  a  com¬ 
plete  representation  can  be  presented.  However,  discussion 
of  the  problem  appears  justified  because  results  for  non¬ 
linear  surging  motion  in  an  Irregular  sea  are  now  available. 

The  problem  is  discussed  for  only  a  long-crested 
sea  and  also  for  only  small  heading  angles.  It  is  hoped  that 
the  reasons  for  broaching  can  be  recognized  by  an  investiga¬ 
tion  within  these  specified  limits. 

THE  EXCITING  FORCES  FOR  THE  SURGING,  SWAYING 
AND  YAWING  MOTIONS _ 

It  is  assumed  that  the  Important  nonlinearities  are 
not  connected  with  the  determination  of  the  exciting  forces 
but  with  the  solutions  of  the  equations  of  motion.  In 
particular,  it  is  assumed  that  the  important  nonlinearities 
are  connected  with  the  displacement  of  the  ship  relative  to 
the  wave  pattern.  This  was  demonstrated  in  part  I  for  the 
surging  motion  in  a  following  sea. 
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No  new  method  Is  used  in  determining  the  exciting 
forces.  Certainly,  there  are  difficult  hydrodynamic  problems 
connected  with  an  accurate  determination  of  these  forces. 
However,  the  aim  of  this  report  is  not  to  treat  these  hydro- 
dynamic  problems  but  to  treat  the  nonlinear  equations  of 
motion.  Therefore,  only  the  Froude-Kryloff  hypothesis  is 
used  for  determining  the  exciting  forces.  It  should  be 
mentioned  that  the  computation  of  the  surging  force  in  this 
way  is  supposed  to  be  a  good  approximation,  while  it  is 
thought  that  the  yawing  moment  and  the  swaying  force  computed 
in  this  way  will  be  too  small.  This  supposition  is  supported 
by  the  fact  that,  for  a  two-dimensional  body  in  a  beam  wave, 
the  computed  swaying  force  is  too  small  if  the  Froude-Kryloff 
hypothesis  is  used  for  such  a  computation.  If  the  swaying 
force  or  yawing  moment  is  computed  by  this  simple  method,  it 
is  estimated  that  the  result  should  be  multiplied  by  a  factor 
of  2. 


Although  the  exciting  forces  resulting  from  the 
Froude-Kryloff  hypothesis  are  used  for  the  numerical  results, 
the  foregoing  remarks  about  such  results  should  be  borne  in 
mind. 


Analogous  to  eq.  3,  the  following  equations  are 
used  for  the  exciting  forces  in  a  regular  quartering  sea 
(  see  Fig .  12) : 


surging  force  = 


swaying  force  = 
yawing  moment  = 


2h 

X 


7T  COSf 


7T  sin^ 
2h 

7 r  sin$ 


(59) 
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It  will  be  permissible  to  replace  cos  f  by  1.0  and 
sin  ip  by  ip  for  a  small  range  of  the  heading  angle  ip . 

The  force  functions  of  ( ^) are  real  functions  for  a 
body  symmetrical  about  the  midship  section .  Furthermore, 
the  functions  f ^  (  j-)  and  f  (£)  will  be  Identical.  For  a 

parabolic  sectional  area  curve,  the  following  expressions 
are  valid: 


f  ( h  =f  (h,  =-a 

xxvL'  yy^v  2 


zz' 


aln(lr)~fe)c°s  (%r) 


m 


(60) 


These  expressions  for  the  exciting  forces  are  used 
for  the  following  numerical  computations. 

For  the  case  in  a  long-crested  irregular  sea, 
eqs.  59  are  replaced  by  the  following: 


Fx  =  surging  force  =  7  cos  ip  ^  cos 


cd2x. 


O  -  CD 


(1- 


L  S 


~)t  +  a(co)J 


fxx  da) 
g2  xx  I 


Fy  =  swaying  force  =  V  sln^  J  cos 


cd2x 


o  -  cd  ( 1 


L  g 


t  i  <r>  <*» 

£ 


(eq.  6l  con't  on  next  pg) 
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00 


F  =yawing  moment  = 
z 


VL  sin  ^1  sin 

0‘ 


i 


d‘ xo  -  m  (1  -  a  t+B(a>) 
g  g 


Vr? *■»  * 

V  6  (61) 

It  is  important  to  note  that  the  following  ratio 
depends  only  on  the  heading  angle  in  an  irregular  long-crested 
sea  as  well  as  in  regular  waves. 


swaying  force 
surging  force 


tan  f 


Nondimensional  spectra  are  computed  assuming  a 
Neumann  spectrum  for  the  long-crested  irregular  sea  and 


plotted  versus  oo 


|  in  Fig. 


13. 


EQUATIONS  OF  THE  SWAYING  AND  YAWING  MOTIONS 

For  investigations  of  coursekeeping  or  steering  in 
calm  water,  the  following  linearized  equations  are  commonly 
used: 

-  V  +  K  -  V  *  -  V  '  V  m 

1^/3  -  Mr('  =  N6S 

As  a  criterion  of  course  stability  in  calm  water, 
the  following  condition  is  known: 

Y/3Nr  +  (rax  "  V  V°  (63) 

Now  it  must  be  decided  what  supplements  are  needed 
for  a  realistic  investigation  of  the  behavior  in  a  severe 
Irregular  quartering  sea.  First  the  following  nondimensional 
expressions  for  the  exciting  swaying  force  and  exciting  yawing 
moment  must  be  added  to  eq.  62: 
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to  the  first  equation: 


pLHV2 

2 


00 

—  sin  ip  f  cosr1 

J  L" 


jd  -  cn  ( 1-  t  +  a  (o>) 

"  O 


to  the  second  equation: 


=  fbc  sin  f 


__„>r  sinf^o  -  .  VT/  , 

£L2HV2  ~2  J  L—  g  zz  L 

2  rr  o  S 

Another  supplement  appears  necessary.  In  finding 
eqs.  6 2  in  ref.  5,  a  constant  speed  of  the  ship  is  assumed. 
However,  a  ship  in  waves  always  will  be  accelerated  or 
decelerated  in  the  longitudinal  direction.  The  influence  of 
such  an  acceleration  is  included  in  eqs.  4  in  ref.  5  and  it 
is  evident  that  this  influence  will  be  important  for  the 
discussions  in  this  report.  The  nondimensional  expression 
V' 

Y  is  used  by  Davidson  in  ref.  5  for  the  acceleration  in  a 

longitudinal  direction.  In  part  I  of  this  report,  the 
acceleration  b  equals  the  right-hand  side  of  eq.  23  if  the 
damping  term  is  neglected.  Using  this  result,  the  non- 
dimensional  acceleration  amounts  to: 


0,(1  .  22)  t  +a(o>) 


/  cd4  f2  (t-)  da> 
r  —  xxvI/ 

I  S2 


.ffll  f “  (t)  d“ 


=  it  -  -L  fcos  [f2xo  -  ®  (1"  t  +  o(co)  /pfiV  (r) 

V2  P  2  J  Ls  “V  gs  L 


The  supplemented  equations  now  can  be  written  as: 

-"Vt-p-  V  +  'vV  *  -  V  Ys6  -Vt 


nz  VU  +  -  iy  -  V  - 

CO 

"•*  +  i  *  Jain  [!%  -  .  d-fW 
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It  is  believed  that  these  supplements  are  suffi¬ 
cient  for  a  discussion  of  broaching. 

The  equations  are  pseudolinear  because  the 
coefficients  of  two  terms  in  each  equation  are  non-constant, 
being  dependent  on  time.  It  is  not  possible  to  linearize 
these  equations  and  still  achieve  the  objectives  of  this 
report.  A  complete  solution  of  these  equations  is  not  known 
and  it  is  not  the  intent  of  this  report  to  develop  such  a 
solution. 

AN  ESTIMATE  OP  COURSE  KEEPING 

IN  A  SEVERE  FOLLOWING  OR  QUARTERING  SEA 

An  attempt  could  be  made  to  replace  the  non¬ 
constant  coefficients  in  the  equations  by  harmonic  functions. 
However,  a  study  of  Pig.  13  shows  that  such  an  approximation 
would  lead  to  a  frequency  for  the  exciting  yawing  moment 
different  than  that  for  the  surging  or  swaying  exciting  force. 
Therefore  a  solution  of  the  equation  would  be  very  complicated 
for  such  an  assumption. 

Another  simplification  of  the  problem  will  be  used 
for  a  rough  estimate.  To  estimate  the  influence  of  the  yaw¬ 
ing  exciting  moment,  a  ship  running  with  a  constant  speed 
determined  by  a  Proude  number  of  about  0.44  is  considered. 

Figure  13  shows  that  for  such  a  run  the  average  period  of  the 
yawing  exciting  moment  will  be  the  largest  possible  period. 

The  influence  of  the  yawing  moment  is  supposed  to  be  the 
largest  for  this  case  (because  the  frequency  of  encounter 
then  equals  zero  for  co  yvfe  = 2  .3) .  Then  the  coefficient  of  the 
yawing  exciting  moment  in  eq.  66  is  replaced  by  a  constant 
coefficient: 
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During  the  same  time  interval  the  surging  and  swaying  excit¬ 
ing  forces  are  unknown.  Perhaps  the  probability  for  dif¬ 
ferent  forces  connected  with  the  moment  of  exp.  67  could  be 
found  by  a  complicated  investigation.  This  will  not  be 
attempted,  however,  and  therefore  the  forces  are  neglected 
for  the  case  considered  here.  In  reality,  the  course  keeping 
will  be  improved  in  some  cases  and  deteriorated  in  other 
cases  by  these  forces.  Therefore,  in  some  cases  the  course 
keeping  will  be  worse  than  estimated  when  neglecting  surging 
or  swaying  exciting  forces. 

Steering  is  necessary  to  avoid  an  unstable  run  of 
the  ship.  The  simplest  assumption  is  chosen  for  steering. 
That  is,  the  rudder  angle  6  is  proportional  to  the  difference 
between  the  heading  angle  ^  to  the  waves  and  the  desired 
heading  angle 

6  =  y(f  -  fQ)  (68) 

The  minimum  factor  7  for  steering,  to  avoid  an  unstable  run 
of  the  ship,  will  be  computed  and  used  as  a  yardstick  for  the 
possibilities  of  steering  in  the  cases  discussed. 

Now,  the  first  very  simplified  case  follows: 

-my/3  +{mx-Yr)i  -  Y  ^  P  =  Y&y(f-fQ) 

-  Nri  =  N6 7(  i/-fo)  —J'a^r  a£f J z ( £)  do>  (69) 

"a"  is  chosen  as  2.25  for  the  average  of  the  l/lOO  most 
severe  cases.  The  spectrum  is  taken  as  in  Fig.  13.  The 
ship's  parameters  are  taken  from  ref.  6  and  ships  A  and  C 
mentioned  in  this  reference  are  chosen.  Then  the  equations 
deliver  the  following  minimum  factor  7: 
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0 

for  ship  A 

for  ship  C 

for 

u2 

gE 

=  b 

Fr=0.45; 

7  min. 

=  0.51 

0.82 

for 

u2 

gE 

=  b> 

Fr=0.44; 

7  min. 

=  0.42 

O.67 

for 

u2 

SC 

Fr=0.43; 

7  min. 

=  0.36 

0.58 

Another  case  is  discussed  before  any  conclusions 


are  drawn  from  these  results.  A  constant  acceleration  of 
the  ship,  and  therefore  also  a  constant  swaying  exciting 
force  is  possible,  over  a  relatively  long  time  interval,  even 
for  a  low  ship  speed.  That  is  shown  in  Fig.  10  for  a  Froude 
number  of  0.142.  Such  an  accelerated  run  of  a  ship  is  con¬ 
sidered  as  the  next  case.  The  yawing  exciting  moment  is 
unknown  during  such  an  accelerated  run.  This  moment  is 
therefore  neglected.  The  same  remarks  mentioned  in  the  pre¬ 
ceding  case  for  neglecting  the  forces  could  now  be  made  for 
neglecting  the  moments.  For  this  second  very  simplified  case, 
the  following  equations  are  used: 

-iyV^-  V/+  K'Yr^"V  =  (70) 

nz^T^+  nz^  “  Np£-Nr^  =  N67^-^o) 


where 


(71) 


The  computations  are  carried  out  for  the  same  ship's 
parameter  as  in  the  preceding  case.  The  following  results  for 
the  necessary  factor  7  are  obtained: 


l-b  V0-1^  t 

f°rj£=l;  F_.-0.lte,  r 


for  ship  A 
8 . 7 ;  7  min .  1.0 

5.8;  7  min.  0.9 


for  ship  C 
2.4 


2.1 
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The  minimum  factor  7  is  much  larger  for  this  case 
than  for  the  preceding  one.  Therefore  it  is  concluded  that 
the  latter  case  is  more  important  if  broaching  is  considered. 


It  coild  be  demonstrated  that  an  accelerated  run 
lasts  even  longer  for  a  smaller  Froude  number  than  for  a 
larger  one.  However,  there  will  be  a  limit.  If  the  Proude 
number  decreases  under  this  limit,  only  a  short  accelerated 
run  will  be  possible.  The  problem  of  course  keeping  will  be 
the  most  severe  for  a  ship  running  with  the  afore-mentioned 
limiting  Froude  number.  There  are  two  reasons  for  this  con¬ 
clusion:  1)  the  time  Interval  for  the  accelerated  run  is 

the  longest  possible;  2)  the  nondimensional  expressions  for 
v/ 

-y-  or  for  the  exciting  forces  are  the  largest  possible. 


The  minimum  factor  7  necessary  for  this  limiting 
Proude  number  will  be  larger  than  the  above  computed  factor 
for  Pr  =  0.142. 


Already  the  factor  7  for  Fr  =  0.142  appears  to  be 
very  large.  Furthermore,  it  appears  as  a  result  of  this 
estimate  that  the  heading  angle  will  grow  rapidly  during  the 
accelerated  run.  This  growth  will  be  intensified  because  at 
the  beginning  of  the  accelerated  run  the  rudder  angle  will  be 
too  small  due  to  a  relatively  quiet  run  of  the  ship  prior  to 
the  accelerated  run. 


The  accelerated  run  lasts  only  a  limited  time 
interval.  This  time  interval  is  estimated  to  be  too  short 
for  broaching  to  occur.  However,  the  heading  angle  will  be 
larger  at  the  end  of  this  time  interval  than  at  the  beginning. 
Broaching  will  occur  if  a  run  with  a  constant  but  increased 
speed  follows  the  accelerated  run.  Although  then  the  surging 
force  is  needed  to  overcome  the  increased  ship’s  resistance, 
a  swaying  force  also  exists  in  this  second  time  interval. 
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The  following  equations  must  be  used  for  this 
second  time  interval  instead  of  eq.  70s 


-^K'yr>*-V  =  V(^o>  -V  / fr  — f dto 


(72) 


nj  -  Npfi  - 

The  following  mimimum  factor  7  computed  from  these  equations 
amounts  to: 


for  ship  A 

for  ship  C 

u2 

ir 

-b 

Fr=0.48;  7  min. 

0.43 

O.78 

u2 

11 

LO|l— * 

F  =0.47;  7  min. 

0.31 

0.55 

This  minimum  factor  7  appears  to  be  too  large  for  ship  C  and 
therefore  it  is  concluded  that  broaching  cannot  be  prevented 
for  ship  C  in  such  seas . 


CONCLUSIONS 


The  course  keeping  of  a  ship  in  a  severe  following 
or  quartering  sea  is  discussed  only  in  a  rough  fashion.  How* 
ever,  the  following  conclusions  appear  to  be  justified: 

1.  The  surging  and  swaying  exciting  forces  are  more 
Important  for  broaching  than  the  yawing  exciting  moment. 

2.  Broaching  probably  will  be  more  severe  for  a  slow  running 
than  for  a  fast  running  ship,  provided  that  the  slow  running 

ship  still  can  be  accelerated  by  the  sea  to  the  dominant  wave 
speed  ( see  Appendix) . 

3.  Broaching  is  more  probable  for  a  ship  course-unstable 
in  calm  water  than  for  a  ship  course- stable  in  calm  water. 

4.  The  suspected  connection  between  broaching  and  the 
nonlinear  surging  motion  is  confirmed  by  conclusion  2. 

5.  The  criterion  for  an  acceleration  of  the  ship  to  the 
dominant  wave  speed  can  be  used  as  an  important  criterion  for 
broaching.  This  criterion  will  be  more  severe  for  a  slow 
running  ship  than  for  a  faster  running  ship. 


R-929 

-39- 


6.  The  method  developed  in  Part  I  can  be  used  for 
determining  this  criterion. 

The  most  interesting  of  these  conclusions  is  the 
second  one.  To  emphasize  this  conclusion,  it  is  added  that 
a  ship  running  with  the  dominant  wave  speed  in  a  severe 
following  irregular  sea  is  extremely  hard  to  steer  and  steer¬ 
ing  control  may  be  lost  at  times.  The  tendency  to  broach 
will  always  be  present,  but  a  severe  yawing  motion  and  tendency 
for  broaching  are  not  necessarily  the  same.  Broaching  is  a 
more  spontaneous  motion  and  connected  with  an  Impetuous 
surging  motion.  Therefore  it  appears  possible  that  the  most 
severe  broaching  tendency  will  occur  for  a  ship  running  with 
a  speed  slower  than  the  dominant  wave  speed. 

RECOMMENDATIONS 

1.  Model  tests  to  validate  the  results  predicted  in  Part  I. 

2.  Determine  the  effect  of  hull  form  on  the  parameters  used 
in  Part  I. 

3.  Determine  the  upper  limit  curve  for  broaching  ( see  Appendix). 

4.  Determine  the  effect  of  hull  form  on  broaching. 

5.  Investigate  the  rolling  motion  connected  with  broaching. 
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NOMENCLATURE 


a 

b 

b 

C 

p  _  _X_ 

Pr-  -  Vg L 
^xx*  Fyy*  ^zz 

fxx(j) ,  \ 

fyy(£) *  I 

f  ,\s  ' 

zz[j) 

g 

h(x,t) 


h 

H 

I 


MxxJ  Myy'  Jzz 
mx  =  m  +  Mxx; 


LaH 


"V  -!Lliz 
2  l2h 


m 

N 


coefficient  for  the  average  i  highest  values 
acceleration 

amplitude  of  acceleration  in  regular  waves 
coefficient  for  the  Neumann  spectrum 
Proude  number 
exciting  forces 


nondimensional  coefficients  for  exciting 
forces  in  regular  waves 

gravitational  acceleration 

displacement  of  the  water  surface  in  the 
vertical  direction 

wave  amplitude  of  a  regular  wave 

draft 

moment  of  inertia 
length  of  the  ship 

added  masses  and  added  moments  of  inertia 


nondimensional  total  masses 


mass 

coefficient  of  damping  force  for  the  surging 
motion 
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V  V  N6 


nondime ns lonal  hydrodynamic  moment  coeffi¬ 
cients 


nz  =  \z  +  *zz 


L4H 


r(o>) 

R 

R 


T 

t,  T,  T 

u 

V 


x,  y,  z 


<JD 


P 

X 

V 

* 

ill 

Yo 

6 

7 

a 


nondimens ional  moment  of  inertia 

spectrum  of  a  long-crested  irregular  sea 
ship  resistance 

a  constant  value  which  corresponds  to  the 
effective  resistance  of  a  ship  driven  by  a 
propeller  running  at  constant  speed 

propeller  thrust 

time  or  time  intervals 

wind  speed 

mean  velocity  for  which  calm  water 
resistance  =  propeller  thrust 

coordinates 

displacement  of  the  surging  motion 

nondimensional  hydrodynamic  lateral  force 
coefficients 

circle  frequency  referring  to  a  fixed 
coordinate  system 

freqency  of  encounter 

density  of  water 

wave  length  of  a  regular  wave 

displacement 

heading  angle  from  the  direction  of  the 
regular  or  long-crested  irregular  sea 

desired  heading  angle 

rudder  angle 

rudder  angle  coefficient  for  a  simple  auto' 
matic  steering 

arbitrary  phase  angle 
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APPENDIX 


In  the  conclusions  of  the  second  section  it  is 
established  that:  "Broaching  will  be  more  severe  for  a  slow 
running  ship  than  for  a  fast  running  ship,  provided  that  the 
slow  running  ship  can  be  accelerated  by  the  sea  to  the 
dominant  wave  speed." 

Figure  14  is  prepared  to  determine  the  limit  for 
an  acceleration  to  a  dominant  wave  speed.  The  curves  plotted 
in  Fig.  14  are  computed  as  described  in  Part  1  for  the  follow¬ 
ing  increment  of  velocity: 

0.40  -  -L-  (73) 

VgL 

0>*L 

=  3*1  can  be  computed  for  this  increment  of  velocity  from 
eq.  53*  Figure  11  then  shows  that  a  run  with  an  Increased 
speed  is  probable  for  an  adequate  time  interval.  Therefore, 
an  acceleration  of  the  ship — and  after  that  a  run  with  an  in¬ 
creased  speed — will  be  possible  if  the  windspeed  is  larger 
than  plotted  in  Fig.  14, 

Figure  14  shows  only  the  limit  for  an  acceleration 
to  a  dominant  wave  speed.  But  it  is  concluded  from  the  dis¬ 
cussions  in  Part  2  that  the  limit  for  an  acceleration  to  a 
dominant  wave  speed  may  also  be  the  limit  for  broaching. 

Figure  14  also  shows  the  acceleration  b  for  the  plotted 
wind  speed. 

From  Fig.  14  it  may  be  concluded  that: 

1.  If  a  ship  runs  with  a  velocity  less  than  the 
limit  speed,  it  cannot  be  accelerated  to  a  dom¬ 
inant  wave  speed  and  broaching  cannot  occur. 

2.  If  a  ship  runs  with  a  velocity  near  to  the  limit 
speed,  broaching  will  not  occur  often,  but  when 
it  does,  it  will  be  very  severe. 

3.  If  a  ship  runs  with  a  velocity  greater  than  the 
limit  speed,  broaching  will  occur  more  frequent¬ 
ly  but  not  as  severely. 
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4.  If  a  ship  runs  with  a  velocity  much  greater 
than  the  limit  speed,  it  will  frequently  be 
accelerated  to  the  dominant  wave  speed,  but 
broaching  may  not  occur  because  the  incre¬ 
ment  of  speed  can  only  be  small  and  because 
the  hydrodynamic  forces  on  the  body  and  the 
rudder  will  be  sufficient  to  maintain  a 
steady  course.  There  must  exist,  therefore, 
an  upper  velocity  limit  for  broaching. 

5.  Broaching  will  not  arise  in  a  mild  sea. 

Therefore  the  limit  curve  for  a  mild  sea  may 
not  be  meaningful. 

The  treatment  in  Part  2  is  not  sufficient  for  deter¬ 
mining  the  limits  mentioned  in  points  4  and  5  above.  A  further 
study  is  therefore  necessary  to  determine  the  limits,  the 
influence  of  the  ship  form  on  these  limits  and  the  effects 
discussed  in  this  paper,  and,  further,  to  confirm  these 
limits  by  model  tests. 
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FIGURE  I.  FUNCTION  FOR  THE  SURGING  EXCITING  FORCE 
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FIGURE  3.  INCREMENT  OF  VELOCITY  FOR  A  HARMONIC  MOTION 
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FIGURE  6ft7.  INCREMENT  OF  VELOCITY  IN  AN  IRREGULAR  FOLLOWING  SEA 
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—  nonlinearity  included 
- NONLINEARITY  EXCLUOEO 
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FIGURE  8. INCREMENT  OF  VELOCITY  IN  AN  IRREGULAR  FOLLOWING  SEA 
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FIGURE  K).  INCREMENT  OF  VELOCITY  IN  AN  IRREGULAR  FOLLOWING  SEA  WITH  DAMPING 
INCLUDED 
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DIRECTION  OF  THE 
LONG  CRESTED  SEA 


FIGURE  12.  LOCATION  OF  SHIP  IN  WAVE 
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FIGURE  14.  LOWEST  WIND  SPEED  AND  LOWEST  ACCELERATION 
NECESSARY  FOR  AN  ACCELERATION  TO  A  SPEED 
OF  Fr=0.40 
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